The Universality of Turbulence in Galactic Molecular Clouds 



Mark H. Heyer^, Christopher M. Brunt ^ 



ABSTRACT 

The universahty of interstellar turbulence is examined from observed struc- 
ture functions of 27 giant molecular clouds and Monte Carlo modeling. We show 
that the structure functions, 5v = foH, derived from wide field imaging of ^^CO 
J=l-0 emission from individual clouds are described by a narrow range in the 
scaling exponent, 7, and the scaling coefficient, Vo- The similarity of turbulent 
structure functions emphasizes the universaility of turbulence in the molecular 
interstellar medium and accounts for the cloud-to-cloud size-line width relation- 
ship initially identified by Larson (1981). The degree of turbulence universality is 
quantified by Monte Carlo simulations that reproduce the mean squared velocity 
residuals of the observed cloud-to-cloud relationship. Upper limits to the varia- 
tion of the scaling amplitudes and exponents for molecular clouds are ~10-20%. 
The measured invariance of turbulence for molecular clouds with vastly different 
sizes, environments, and star formation activity suggests a common formation 
mechanism such as converging turbulent fiows within the diffuse ISM and a lim- 
ited contribution of energy from sources within the cloud with respect to large 
scale driving mechanisms. 

Subject headings: ISM: clouds - ISM: kinematics and dynamics -ISM: molecules 
- ISM: structure 



1. Introduction 



Turbulent, non-laminar gas fiows are a ubiquitous feature within all phases of the inter- 
stellar medium (ISM). Therefore, accurate descriptions of interstellar turbulence are essential 
to meaningful understanding of ISM dynamics and star formation. An important statistical 
description of fiuid dynamics is the generalized velocity structure function, 

S,il)=<\v{r)-v{r + l)\^> (1), 
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where I is the spatial displacement between two cells within a 3 dimensional volume, p is the 
order of the structure function, and the averages arc taken over the volume of the fluid. Over a 
finite spatial range, the structure functions may be described as a power law, Sp{l) ~ l'^^^^ but 
is often re-framed as an equivalent linear expression by taking the pth root, Sp^iy^^ = 6v = 
Vol'^, where 7 = ({p) /p is the scaling exponent and Vo is the scaling coefficient. The structure 
function provides a concise description of the spatial coherence of velocity differences within 
a given region. Such differences can arise from both systematic motions (rotation, collapse, 
outflows) and random motions due to turbulent gas flows. Within interstellar clouds, most 
velocity differences are due to turbulence. 

One of the most cited and influential studies on interstellar turbulence is Larson (1981) 
that identified a power law relationship between the global velocity dispersion, Av (km s~^), 
and cloud size, L (pc), of molecular clouds from values taken from the literature, Av — CLF 
where C — 1.1 and F = 0.38. Using a more homogeneous set of cloud data from the 
Massachusetts-Stony Brook Galactic Plane Survey, Solomon et al. (1987) found a similar 
correlation with a comparable scaling coefficient {C = 1.0) but a steeper scaling exponent 
(r = 0.5). These relationships rely on global velocity dispersions that are distinguished 
from velocity differences codified in the structure function of equation (1). As the sample 
clouds are distributed throughout the Galaxy, these do not comprise a singular fluid volume. 
Therefore, the connection of these cloud-to-cloud size-line width relationships to velocity 
structure functions of individual clouds would seem remote unless one assumes that all 
clouds have approximately the same values for Vo and 7. Larson (1981) showed examples 
of a similar scaling law within clouds by using sizes and velocity dispersions derived from 
molecular tracers with different excitation requirements (see also Fuller & Myers 1992). 
However, the dense regions comprise a small fraction of the mass and volume of a molecular 
cloud so the measured velocity dispersions within the spatial scales of CS or NH3 emissions 
may not be representative of velocity differences over comparable scales but within a lower 
density substrate. Size-line width relationships derived from clump identiflcation algorithms 
for individual clouds exhibit a large degree of scatter with a range of values for 7 (Carr 
1987; Stiitzki & Giisten 1990; Falgarone, Puget, & Pereault 1992) or no relationship at 
all (Loren 1989; Simon etal 2001). The absence of a consistent correlation between size 
and line width of embedded structures within a given cloud is due to the limited dynamic 
range of sizes that can be identified by such algorithms and the superposition of emission 
from disconnected features along the line of sight (Ballesteros-Paredes & Mac Low 2002). 
Using structure functions of velocity centroid images, Meisch & Bally (1994) showed that 
the scaling exponents are similar for a sample of 4 molecular clouds but did not consider the 
variation of the scahng coefficient. Brunt (2003) provides the most convincing evidence for 
similar scaling laws within molecular clouds using Principal Component Analysis (PCA) as 
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a tool to recover the true structure function for a given cloud (Brunt & Heyer 2002). For 
each ^^CO or ^^CO spectroscopic data cube of a molecular cloud, a set of 5v,l points are 
determined from the eigenvectors and eigenimages. When the PCA measurements for all 
clouds are combined onto a single plot, these define a nearly co-linear set of points. Such 
a correlation necessarily results from narrow distributions of the scaling exponent, 7, and 
coefficient, Vo, for this sample of clouds. 

In this Letter, we extend the study of Brunt (2003) to demonstrate that Larson's cloud- 
to- cloud scaling law is explained only if the structure functions for individual clouds are 
nearly identical. PCA-based 6v, I relationships arc presented that demonstrate the same 
functional form for structure functions for molecular clouds that span a wide range in size 
and environmental conditions. Monte Carlo models are constructed that place upper limits 
to the variation of the scaling coefficient and exponent. Finally, we discuss the consequences 
of an invariant turbulent spectrum in the context of the formation of interstellar molecular 
clouds, sources of turbulent energy, and star formation. 



2. The Composite Structure Function 

Following Brunt & Heyer (2002), PCA is apphed to spectroscopic data cubes of ^^CO 
J=l-0 emission from molecular clouds that are part of recent wide field imaging surveys 
at the Five College Radio Astronomy Observatory (Heyer et al. 1998; Brunt & Mac Low 
2004) or targeted studies of individual giant molecular clouds. Heyer & Schloerb (1997) 
and Brunt (2003) show there is little difference in the Sv, I relationships derived from ^^CO 
emission and the lower opacity ^^CO emission. For each cloud, a power law is fit to the 
Sv,l points to determine the PCA scaling exponent, apcA, and coefficient, Vo- For the 
sample of 27 molecular clouds, the mean and standard deviation for the scaling exponent 
are 0.62 and 0.09 respectively. Based on models with fittle or zero intermittency, this PCA 
scaling exponent corresponds to a structure function exponent equal to 0.49 ± 0.15 (Brunt 
et al. 2003). The mean and standard deviation of the scahng coefficient are 0.90 km s~^ 
and 0.19 km s~^. These rather narrow distributions of 7 and Vo re-emphasize the results of 
Brunt (2003) that there is not much variation in the structure function parameters between 
molecular clouds. In Figure 1, we overlay the PCA 6v, I points from the sample of clouds. 
The composite points reveal a near-identical form of the inferred structure functions. The 
solid line shows the power law bisector fit to all points, 6v = (0.87 ± 0.02)/°-65±°-°^ This 
PCA derived exponent corresponds to a structure function scaling exponent of 0.56±0.02. 

The global velocity dispersion of each cloud and the cloud size are determined from 
the scales of the first eigenvector and eigenimage respectively. Basically, the global velocity 
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dispersion, Av, is the value of the velocity structure function measured at the size scale, 
L, of the cloud. These points, marked as filled circles within Figure 1, are equivalent to 
the global values used by Larson (1981) and Solomon et al. (1987) that define the cloud- 
to-cloud size-line width relationship. A power law bisector to this subset of points is Av = 
(0.96 ± 0.17)L'''^^^°''^^. The similarity of this cloud-to-cloud relationship with that of the 
composite points is a consequence of the uniformity of the individual structure functions. 
Within the quoted errors, it is also similar to the cloud-to-cloud size-hne width relationships 
- 7 ~ r and Vo ~ C. Therefore, Larson's global velocity dispersion versus cloud size scaling 
law follows directly from the near identical functional form of velocity structure functions for 
all clouds. If there were significant differences of 7 or Vo between clouds, then the cloud- 
to-cloud size-line width relationship would exhibit much larger scatter than is measured by 
Larson (1981) and Solomon et al. (1987). 



3. The Degree of Turbulence Universality 

The cloud-to-cloud size-line width relationships measured by Larson (1981) and Solomon 
et al. (1987) and the composite structure functions shown in Figure 1 do exhibit some degree 
of scatter about the fitted lines. The scatter is quantified by the mean square of the velocity 
residuals, cr^^^, for each data set where 

N is the number of clouds in the sample, C and F are the parameters derived by fitting 
a power law to the observed A^;, L points. The value for a^^g for the sample of clouds in 
Larson (1981) using only the -"^^CO and ^^CO measurements is 1.41 fcm^s"^. The Solomon 
et al. (1987) sample is a larger, more homogeneous set of clouds and therefore, provides a 
more accurate measure of the variance within the cloud-to-cloud size-line width relationship. 
The corresponding cr^j,^ is 0.88 /cm^s"^. The value of a'^f^^ for the Av, L points in Figure 1 is 
1.93 krn^s"^ and 0.35 krn?s~'^ for the composite collection of 5v, I points. 

The measured scatter, described by 0"^^^, of the size-line width relationships is a critical 
constraint to the degree of invariance of turbulence within the molecular interstellar medium. 
The scatter arises from several sources. There are basic measurement errors in the global 
velocity dispersion owing to the velocity resolution of the measurements and the cumulative 
statistical error of the individual spectra. Deriving cloud sizes from complex projected 
distributions of the molecular gas may also introduce some scatter. These measurement 
errors are rarely shown in any cloud size-line width plots. A secondary source of scatter is 
limited or biased mapping of the molecular cloud. If a given map was limited in angular 
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extent and centered on a region within the cloud that is actively forming stars then the 
measured "global" velocity dispersion may be broader due to localized, expanding motions 
from HII regions or protostellar outflows. Such enhanced velocity dispersions from biased 
regions may not represent the global velocity dispersion of the entire giant molecular cloud, 
that is typically quite extended with respect to star formation sites within the cloud. An 
additional source of the observed scatter is the uncertainty of the distance to each cloud in 
the sample. Distances to molecular clouds, and correspondingly, cloud sizes, are generally 
not known to precisions better than 25%. Finally, and most important for the subject of 
this study, true differences in the scaling coefficient and exponent would also contribute to 
the observed scatter in the size-line width relationship. It is this component that we wish to 
constrain as it defines the degree to which turbulence is invariant in the ISM. 

To gauge the universality of turbulence within the molecular interstellar medium, a 
simple, Monte Carlo model is constructed to place limits on the variance of the scahng 
coefficients and scaling exponents of individual clouds. The structure function for each 
model cloud is defined by the parameters, t'o =< t'o > +ey and 7 =< 7 > where e„ 
and are drawn from gaussian probability distributions with standard deviations of 
and (7^ respectively. Both (j^, and implicitly contain contributions from the measurement 
errors, biased mapping, and true variations in turbulence. In the limit of absolute universal 
turbulence, infinite precision, and unbiased imaging, ay = a^=0. A size, Lo, is assigned to 
each cloud from a uniform probability distribution such that < log{Lo) < 2, corresponding 
to cloud sizes between 1 and 100 pc. A global velocity dispersion, Av, is determined by 
evaluating the structure function at the assigned size of the cloud, 

Av = (< > +e,)L<^>+^- (3) 

Once Av is determined from equation (3), a random component, ei, is added to the cloud 
size, L = Lo(l + e^) to replicate an uncertainty in cloud distances where e^^ is drawn from a 
uniform probability distribution within the percentage range, ±(Td/£). Simulations are run 
for ao/D—O, 0.1, 0.25, and 0.5. The Solomon et al. (1987) data set is used as the primary 
observational constraint so N=272. Following the results in Figure 1, < Vo >=0.9 km s~^, 
and < 7 >=0.5. For each set of (T„ and ay parameters, is calculated, 

a{ay,ay)^ km s (4), 

where Ci and Fi are determined from a bisector fit to the Av, L points for a single realization. 
To reduce the statistical errors of the simulation, we calculate the mean value of cr^ ((7^,(7^) 
from 500 realizations. 

The Monte Carlo results are shown in Figure 2 where (7„ and are normalized hy < Vo > 
and < 7 > respectively to reflect the fractional variation. The heavy solid hne shows the 
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value of kw?s~'^ from the Solomon et al. (1987) sample. This contour defines the 

locus of a^/ < Vo >, cr^/ < 7 > points that reproduce the observed scatter in the Solomon 
et al. (1987) size-line width relationship. In the unlikely extreme limit that a^=0, then the 
most (T„ can vary is 18-23% about the observed value of 0.9 km s^^. Conversely, if (Tt,=0, then 
can vary by less than 9-12% about its fiducial value of 0.5. More realistically, cr^ ^ and 
^ 0, so the percent variations for both parameters are ~8-12%. We emphasize that these 
are upper limits to the true variations between clouds as (7„ and also include measurement 
errors and biased mapping effects that are likely present in all cloud-to-cloud size-line width 
relationships. 



4. Implications of Invariant Interstellar Turbulence 

The upper limits to variations of the structure function parameters are quite small when 
one considers the large range of molecular cloud environments. With few exceptions, giant 
molecular clouds are sites of massive star formation. The massive stars can, in principle, 
affect gas dynamics over the extent of a cloud by enhancing the UV radiation field, driving 
HII region expansions and stellar winds, and are the progenitors of supernova explosions. 
The star formation activity in smaller clouds is generally limited to the birth of low mass stars 
whose impact on cloud dynamics by protostcllar winds is highly localized and small with 
respect to the integrated energy input from massive stars. For the sample of molecular clouds 
used in this study, the far infrared IRAS point source luminosities range from 20 Lq (B18, 
Heiles' Cloud 2, L1228) to 7.6 xlO^ Lq (NGC 7538). Despite the large differences in the 
amount of internal energy injection from newborn stars, the functional forms of the turbulent 
structure functions for all clouds are similar. Either all clouds coincidentally redistribute this 
internal energy into turbulent, random motions described by structure functions with the 
same slopes and amplitudes via some self-regulatory process or these internal energy sources 
are small compared to an external energy reservoir that is common to all regions in the 
Galaxy. Brunt ct al. (2004) show that most of the turbulent energy of a molecular cloud 
resides on the largest scales and that the cloud dynamics are more readily accounted by large 
scale driving of turbulence. The universality of turbulence described in this study provides 
additional evidence for large scale driving sources within the molecular interstellar medium. 
The narrow range of turbulent fiow parameters may also refiect the necessary conditions 
that facilitate the development of molecular clouds. Regions with extreme, high values of Vo 
may be highly overpressured with respect to self-gravity and expand to larger, more diffuse 
configurations with less effective self-shielding to sustain significant molecular abundances. 

Larson (1981) presaged the results from recent numerical simulations that molecular 
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clouds are short lived (<10^ years), transient objects (Ballesteros-Paredes, Hartmann, & 
Vazquez-Semadeni 1999) that must be continually reassembled from the diffuse gas compo- 
nent. His suggestion that molecular clouds result from thermal instabilities within shocks of 
colliding atomic gas streams is similar to recent numerical simulations of compressible tur- 
bulence that show the development of high density regions directly from the shocks (Hunter 
et al. 1986; Elmegreen 1993; Ballesteros-Paredes, Vazquez-Semadeni, & Scalo 1999). Such a 
common dynamical origin of molecular clouds in the Galaxy may account for the measured 
near invariance of turbulence. 

The results presented in this study apply to the low density gas substrate that comprises 
most of the mass and volume of a molecular cloud. Turbulent properties may indeed differ 
between clouds within the high density, localized, supercritical core regions. In distributed, 
low mass star forming regions, the non-thermal motions within the dense gas are subsonic 
(Benson & Myers 1989). Within the massive cores typical of clustered star forming regions, 
the observed velocities are supersonic (Pirogov et al. 2003). The identification of processes 
responsible for such divergent, dense gas configurations remains one of the primary challenges 
to descriptions of star formation (Klessen, Heitsch, & Mac Low 2000; Padoan & Nordlund 
2002; Vazguez-Semadeni, Ballesteros-Paredes, & Klessen 2003; Shu et al. 2004; Mac Low & 
Klessen 2004). 

5. Summary 

We have examined the degree to which interstellar turbulence is universal within the 
molecular gas component of the Galaxy by comparing the measured structure functions for 
27 giant molecular clouds. Despite the large differences in cloud environments and local 
star formation activity, the structure functions are described by very narrow ranges of the 
scaling exponent and scaling coefficient. The degree of universality is further constrained 
by Monte Garlo simulations that replicate the observed scatter in the Larson scaling law 
that describes the relationship between global velocity dispersion and cloud size. The near 
invariance of turbulence within the molecular interstellar medium suggests a common for- 
mation mechanism of molecular clouds such as shocks due to colliding gas streams of diffuse, 
atomic material as originally suggested by Larson (1981). It also implies that the energy 
contribution from internal sources such as stellar winds and expanding HII regions may be 
small with respect to a common, external component. 

We acknowledge useful comments from Richard Larson, John Scalo, Bruce Elmegreen, 
Paolo Padoan, and Enrique Vazquez-Semadeni. This work was supported by NSF grant 



-8- 



AST 02-28993 to the Five College Radio Astronomy Observatory. 

REFERENCES 

Ballesteros-Paredes, J.E., Vazquez-Semadeni, E., & Scalo, J. 1999, ApJ, 514, 286 
Ballesteros-Paredes, J.E., Hartmann, L., & Vazquez-Semadeni, E. 1999, ApJ, 527, 285 
Ballesteros-Paredes, J.E. & Mac Low, M. 2002, ApJ, 570, 734 
Benson, P.J. & Myers, P.C. 1989, ApJS, 71, 89 
Brunt, CM. & Heyer, M.H. 2002, ApJ, 566, 276 
Brunt, CM. 2003, ApJ, 584, 293 

Brunt, CM., Heyer, M.H., Vazquez-Semadeni, E., & Pichardo, B. 2003, ApJ, 595, 824 
Brunt, CM. & Mac Low, M. 2004, ApJ, 604, 196 

Brunt, CM., Heyer, M.H., Zivkov, V., & Mac Low, M. 2004, ApJ, in press 
Carr, J.S. 1987, ApJ, 323, 170 
Elmegreen, B.C. 1993, ApJ, 419, L29 

Falgarone, E. Puget, J.-L., Pereault, M. 1992, A&A, 257, 715 
Fuller, G.A. & Myers, P.C 1992, ApJ, 384, 523 
Heyer, M.H. & Schloerb, F.P. 1997, ApJ, 475, 173 

Heyer, M.H., Brunt, C, Snell, R.L., Howe, J.E., Schloerb, F.P., & Carpenter, J.M. 1998, 
ApJS, 115, 241 

Hunter, J.H., Sandford, M.T., Whitaker, R.W., & Klein, R.I. 1986, ApJ, 305, 309 
Larson, R.B. 1981, MNRAS, 194, 809 
Loren, R.B., 1989, ApJ, 338, 925 

Klessen, R.S., Heitsch, F., & Mac Low, M. 2000, ApJ, 535, 887 

Mac Low, M., & Klessen, R.S. 2004, Reviews of Modern Physics, 76, 125 



-9- 

Padoan, P. & Nordlund, A. 2002, ApJ, 576, 870 

Pirogov, L., Zinchenko, I., Caselli, P., Johansson, L.E.B., & Myers, P.C. 2003, A&A, 405, 
639 

Shu, F.H., Li, Z-Y, & Allen, A. 2004, ApJ, 601, 930 

Simon, R., Jackson, J.M., Clemens, D.P., Bania, T.M., & Heyer, M.H. 2001, ApJ, 551, 747 
Solomon, P.M., Rivolo, A.R., Barrett, J., & Yahil, A. 1987, ApJ, 319, 730 
Stiitzki, J. & Glisten, R. 1990, ApJ, 356, 513 

Vazquez, E., Ballesteros-Paredes, J.E. & Klessen, R.S. 2003, ApJ, 585, L131 



This preprint was prepared with the A AS IM^jX macros v5.2. 



-10- 



T 1 1 1 — I I I I I 1 1 1 1 — I I I I I 1 1 1 1 — I I I I I 1 1 1 1 — I I I I 




I (pc) 

Fig. 1. — The composite 5f , / relationship from PCA decompositions of ^^CO J=l-0 imaging 
observations of 27 individual molecular clouds. The small scatter of points attest to the near 
invariance of interstellar turbulence within molecular clouds that exhibit a large range in 
size, environment, and star formation activity. The large filled circles are the global velocity 
dispersion and size for each cloud derived from the first principal component. These are 
equivalent to the global velocity dispersion and size of the cloud as would be measured in 
the cloud-to-cloud size-line width relationship (Larson 1981; Solomon et al. 1987). The hght 
solid line show the bisector fit to all points from all clouds. The heavy solid line shows the 
bisector fit to the filled circles exclusively. The similarity of these two power laws explains 
the connection of Larson's cloud-to-cloud scaling law to the structure functions of individual 
clouds. 
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Fig. 2. — Images of from the size-line width fits for an ensemble of clouds with varying 
turbulent properties parameterized by (y^l< Vo > and o'^/< 7 > for different distances 
uncertainties (o"d/-D=0, 0.1,0.25, and 0.5). The light contour values of o"^ are 0.5, 1.0, 1.5 2, 
2.5, 3.0 kw?s~^. The heavy contour in each plot shows the loci that correspond to (Tpj,^=0.88 
km^ s~^ from Solomon et al. (1987) and provides an observational upper limit to cloud-to- 
cloud variations of Vo and 7. These loci demonstrate that the structure function parameters, 
7 and Vo-i can not vary by more than ~10-20%. 



